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Abstract 

In this paper, we will give convergence analysis for a family of 14-node elements 
which was proposed by I. M. Smith and D. J. Kidger in 1992. The 14 DOFs are 
taken as the value at the eight vertices and six face-centroids. For second-order 
elliptic problem, we will show that among all the Smith-Kidger 14-node elements, 
Type 1, Type 2 and type 5 elements can get the optimal convergence order and 
Type 6 get lower convergence order. Motivated by the proof, we also present a 
new 14-node nonconforming element. If we change the DOFs into the value at the 
eight vertices and the integration value of six faces, we show that Type 1, Type 
2 and Type 5 keep the optimal convergence order and Type 6 element improve 
one order accuracy which means that it also get optimal convergence order. 

Keywords: Nonconforming element; Brick element; 14-node element; Second- 
order elliptic problem; Smith-Kidger element 



1 Introduction 

Until now, there are many three-dimensional brick elements. For conforming case, 
the trilinear element, 27-node element and Seredipity elements are well known. For 
nonconforming case, Rannacher-Turek ^ presented the rotated Qi nonconforming ele- 
ment, Douglas-Santos-Sheen- Ye [2j introduced the nonconforming finite element using 
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mental Research Funds for the Central Universities" . Also partially supported by NRF of Korea (Nos. 
2011-0000344, FOl-2009-000-10122-0). 
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only six values at the centroids of the faces as degrees of freedom, and Park- Sheen pre- 
sented a Pi-nonconforming finite element on cube meshes which has the lowest degrees 
of freedom [6j . Wilson also defined a linear-order nonconforming brick element [H [15] 
with 11 DOFs whose basis consists of trilinear polynomials plus {1 — x"^ , 1 — y'^ , 1 — z"^} 
on K = [—1, 1]^. This element obtains an 0{h) convergence rate in energy norm, (see 
m Page 217, Remark 4.2.3]) 

To get a higher order of convergence for error estimates. Smith and Kidger [12] 
successfully developed three-dimensional brick elements of 14 DOFs. They investigated 
six most possible 14 DOFs elements systematically considering the Pascal pyramid, 
and concluded that the type 1 (as well as type 2) and type 6 elements are success- 
ful ones. The type 1 element add the span of four nonsymmetric cubic polynomials 
{xyz, x'^y,y'^z, z'^x} while the type 6 element the span of {xyz,xy'^z'^,x'^yz'^,x^y'^z} to 
P2- Only recently a new nonconforming brick element of fourteen DOFs with quadratic 
and cubic convergence in the energy and L2 norms is introduced by Meng, Sheen, Luo, 
and Kim [5], which has the same type of DOFs but has only cubic polynomials added 
to P2- And then, the authors compared these 14-node elements numerically, see |1]. 
Numerical tests showed that at least for second-order elliptic problems Meng-Sheen- 
Luo-Kim and part of Smith-Kidger elements are convergent with optimal order or with 
lower order. 

Convergence analysis for Meng-Sheen-Luo-Kim element was done in |5j and is fairly 
easy because it satisfies the patch test of Irons [3J, which implies that a successful Pk- 
nonconforming element needs to satisfy that on each interface the jump of adjacent 
polynomials be orthogonal to Pk-i polynomials on the interface. Unfortunately, it was 
found in mathematics that the patch test is neither necessary nor sufficient, see [9] 
and the reference therein. As showed in this paper, the Smith-Kidge elements can 
only pass lower order patch test or can not pass it, but give optimal convergence order 
from our numerical results or lower convergence order. Thus, the convergence analysis 
for Smith-Kidger element is quite different and complex. For the convergence of the 
nonconforming element which can not pass the patch test, Stummel, Zhongci Shi etc. 
have done a lot of work, see [El [H [101 E] • 

In this paper, we will give convergence analysis for Smith-Kidger element for the 
second-order elliptic problem. We show that although the patch test can not be satisfied. 
Type 1, 2 and 5 Smith-Kidger element can also get optimal convergence order, and Type 
6 element loses one order of accuracy. Furthermore, we also present a new brick element 
with the same DOFs, which is also convergent with optimal order. Finally, if the value 
at the eight vertices and the integration value of six faces are taken as the DOFs, then 
we can show that Type 1, 2, 5, 6 elements and the proposed new element can get 
optimal convergence order, which implies that Type 6 element improves one order of 
accuracy. 

The paper is organized as follows. In section 2, we will introduce Smith-Kidger 
element and give the basis functions firstly. In section 3, we define the interpolation 
operator and present convergence analysis for Type 1 Smith-Kidger element. In section 
4, we will analyze the other elements and present the corresponding error estimates 
very briefiy. In section 5, a new 14-node brick element is proposed. Finally, in section 
6, we will conclude our results. 
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2 The quadratic nonconforming brick elements 

Let K = [—1,1]'^ and denote the vertices and face-centroids by X^-, 1 < j < 8, and 
1 < A; < 6, respectively, (see Fig. [1]) 



Vy 


^ 1\/T 






Me' 
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• Ml 






'Mi 





Figure 1: Vj denotes the vertices, i = 1,2,..., 8, and denotes the face-centroid, 
k = 1,2,. ..,6. 

Smith and Kidger [12] defined the following six 14-node elements: 





= P2 


:k)( 


33 Span{a??^, x^y, y^z, z^x}, 


(la) 


p{2) 


= P2 


:k)( 


33 Spaii{xyz, xjf, yz^, zx"^}, 


(lb) 


p{3) 


= P2 


:k)( 


33 Span{a?2^, x^,y^, z^}, 


(Ic) 


p{4) 


= P2 


:k)( 


3) Spaii{xyz, x^yz, xy^z, xyz^}, 


(Id) 


p{5) 


= P2 


:k)( 


33 Span{xyz, x^y + xy^ , y^^ + yz^ , 'z^x + ^x^}. 


(le) 


p{6) 


= P2 


:k)( 


3) Span{a?y2, xy^z^ ,x^yz^ , x^y^z^}, 


(If) 



whose DOFs are the function values at the eight vertices and the six face-centroids. 
They reported that Type 3 element fails and is inadmissible. We also remark that 
Type 4 element is also inadmissible since (x^ — l)yz G vanishes at all these points. 
In we observe that Type 1 (and 2) and Type 5 elements give optimal convergence 
results both in and norms at least for the second-order elliptic problems, while 
Type 6 element loses one order of accuracy in each norm. 

In what follows, we will give the error estimate for the first type Smith-Kidger 
element detailedly and the error estimates for other types can be done completely 
similarly and are stated very briefiy. 

Firstly, the basis function of first type 

01 = -— + (x'^ + lf + z^) + 77 (xy + xz + yz + xyz + x^y + y'^z + z^x) , 

lo lo 8 

02 = -T^ + (^^ + + + o {xy -xz -yz -xyz + x^y -y^z + z^x) , 
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(P9 = ^ + + -[x -y - z ) - -Z X, 
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—xy + xz — yz — xyz — x^y + y^z + z^x) , 
—xy — xz — yz + xyz — x^y — y^z-\- z^x) , 
—xy — xz + yz — xyz + x^y + y^z — z^x) , 
—xy-\- xz — yz + xyz + x^y — y^z — z^x^ , 



xy — xz — yz + xyz — x^y + ifS — 'z^x 



xy + xz + yz — xyz — x^y — y^z — z^x) 
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4 2 

1 1 



014 = I - 4 (^^ - - + 2^^^- 

Assume that G is a parallelepiped domain with boundary F. Let {S^h)h>o be 
a regular family of triangulation of f2 into parallelepipeds Kj, j = 1,2, ... , A^k, where 
h = maxKg^-^ with = diam(K). For each K e let Tk : K ^ be an 
invertible affine mapping such that 

K = rK(K), 

and denote 0k = o Tj^^ : K ^ M for all e P^X' whose collection will be designated 

by 

Pk = Span{0K, e P^'l}. 
Let AV and Np denote the numbers of vertices and faces, respectively. Then set 



yh = {Vi,v2,...,Vj 



Nv 



the set of all vertices of K G 



= {Fi, F2, . . . , Fnp : the set of all faces of K G ^h}, 

= {Ml, M2, . . . , Mnj, : the set of all face-centroids on ^h}- 

We consider the following nonconforming finite element spaces: 

^'rf'' = {(p:n^R\ 01 K e PkVK G ^, is continuous at all Vj G fh, Mk G ^h}, 
= {0 G 0(y) = VFj G n F and 0(M) = VM^ G n F}. 
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3 The interpolation operator and convergence anal- 
ysis 

In this section we will define an interpolation operator and analyze convergence in 
the case of Dirichlet problem. The case of Neumann problem is quite similar and the 
results will be omitted. 

3.1 The second order elliptic problem 

Denote by (-,■) the L'^{n) inner product and (/, v) will be understood as the duality 
pairing between H^^{Q) and HqIQ), which is an extension of the duality paring between 
L^(0). By II ■ II and | ■ |jt we adopt the standard notations for the norm and seminorm 
for the Sobolev spaces H^{Q). Consider then the following Dirichlet problem: 

- V ■ (aVu) + /3u = f, n, (2a) 
u = 0, r, (2b) 

with a = {ajk),ajk,^ G L^{Q),j,k = 1,2,3,0 < a,|^p < ^*Q:(x)| < a*\^\^ < 
oo,^ G M^,/3(a;) > 0, x E fl, and / G H^{fl). We will assume that the coefficients are 
sufficiently smooth and that the elliptic problem (|2]) has an i/^(f2)-regular solution such 
that ll^lls < C||/||i. The weak problem is then given as usual: find u G Hq{Q) such 
that 

aiu,v) = if,v), veH'.in), (3) 

where a : Hq{Q) x Hq{Q) — > M is the bilinear form defined by a{u,v) = {aVu, Vv) + 
{/3u,v) for all u,v E HI{VL). The nonconforming method for Problem ([2]) states as 
follows: find Uh G such that 

cih{uh. Vh) = if, Vh), Vh G ^'^0, (4) 

where 

ah{u,v) = ^ a-K{u,v), 

with Ok being the restriction of a to K. 

Notice that in order to have point values defined properly we need to recall Sobolev 
embedding theorem 

1 TYl 

W'^'PiQ) C\Q), if - - < 0. 

p a 

Thus we should have p > y. For a given cube K G .5^, define the local interpolation 
operator Hk : W^'P(K) D H^{n) — > ¥^1, p > 3, by 

for all vertices Vi and face-centroids Mj of K. The global interpolation operator Uh. 
W^'P{Q) n Hq{^1) — )■ ^"^Q is then defined through the local interpolation operator IIk 
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by H/jIk = Hk for all K G Since 11^ preserves P2 for all K G it follows from 
the Bramble-Hilbert Lemma that 

Ke-^-fc KG,^h (5) 

G nHl{n),i<k< 3. 

We now consider the energy-norm error estimate and first consider the following 
Strang lemma [T3] . 



Lemma 1. Let u G H^iVl) and Uh G jV^^ he the solutions of Eq. ([3D and Eq. (jlj), 
respectively. Then 

II II ^ / • f II II , \ah{u,w) - {f,w)\ \ 

\\u-Uh\\h<c< mf \\u-v\\h+ sup — >. (6) 



Assume sufficient regularity such that u G H^. Due to ([5]), the first term in the 
right side of (E]) is bounded by 

inf \\u - v\\h < \\u -UhuWh < ch''\u\Hs+i(n)A < s <2 (7) 

Now let us bound the second term of the right side of ([6]) which denotes the con- 
sistency error. For a given cube K G =5^, denote by F^~^ and the face of K with 
outward normal (1,0,0) and (—1,0,0), respectively. By the same way, we denote the 
other faces by , F^~ , F^~^, F^~. Thus, integrating by parts elementwise, we have 

ahiu,w) - {f,w) = ^ (i/^aVM,M;)aK 

= 2 - / ly^ctWuwds + 2. / v^aWuwds 

+ / u^cxVuwds 

'■= Ex + Ey + Ez 

where is the transpose of the unit outward normal to K. 
Before proceeding, we need the following lemma. 

Lemma 2. By Fk denote the face containing the centroid Mk and by V^'' ,j = 1,2,3,4, 

denote the vertices on the surface Fk. If p & '^s^ki Pi^f'') = 0,i = 1,2,3,4, and 
p{Mk) = 0, then 

''p{x,y,z)ds = 0, /c = l,2,...,6. (8) 
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Proof. Without loss of generality, we assume that Mi = (1,0,0). In this case, we have 
p e and ±1, ±1) = 0, 0) = 0. We need to prove that 

/ pil,y,z)dydz = 0. (9) 

JFi 

It follows from p{l, ±1, ±1) = that 

p{l, y, z) = h{y, z){y' - 1) + hiy, z){z^ - 1), G PiiM.'), j = 1, 2. (10) 

Set 

lj{y,z) = ajy + bjZ + Cj, j = 1,2. 
Then p{l, 0, 0) = implies that Ci + C2 = 0, which reduces ffTOj) to 

p{l, y, z) = {aiy + hiz){y'' - 1) + (aay + h2z){z^ - 1) + ci{y'' - z"). 

Since 

P5l-U=i = Span{l, z, y'^, yz, z^ , y'^z}, (11) 
and notice that p G P^^' have 

Ol = 02 = &2 = 0, 

which leads to 

p(l, y, z) = hz{y^ - 1) + ci(y2 - ^2). (12) 
By a direct computation using f|T2|) . Eq. ([8]) holds. This completes the proof. □ 



This lemma implies that Type 1 element can pass lower order patch test (test 
functions are in Pq not Pi), which will lead to a convergence solution for the second 
order elliptic problems, but the convergence order is not optimal. To bound Ex,Ey,Ez, 
we also need some interpolation operators. 

3.2 Some interpolation operators 

For the reference element K = [—1, 1] x [—1, 1] x [—1, 1], we consider the interpo- 
lation problem on the face of F^"*". The interpolation points are (1,1,1), (1,1,-1), 
(1, —1, 1), (1, —1,-1), (1, 0, 0) which are the four vertices and the centroid of F^, and 
the interpolation space is Span{l, 'z, yz, z^}. This interpolation problem is poised and 
the basis functions are 

1_ 1^ 1_. 1^2 

1^ 1^ 1_ 1^2 

^2 = -y--z--yz+-z , 

1^ 1^ 1^ 1^2 

^3 = --y + -z--yz + -z , 
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Thus for a function / defined on the face F^, the interpolation polynomial can be 
written as 

f^+f = /(I, 1, 1)^1 + /(I, 1, -1)^2 + /(I, -1, 1)^3 + /(I, -1, -1)^4 + /(I, 0, 0)^5 

For a given Kg we can define the interpolation operator Tp' by Tp' = Tp^ o Tj^^ 
Similarly, we can also define the interpolation operators on the other five faces of K 
and denote them by T^' ,Tp' ,Tlr ,T^+ ,T^- , respectively. Here we remark that the 
interpolation space of Tp~ is same with T^^, the interpolation space of Tp~ and Tj'r is 
Span{l, X, z, xz, x^}, and the interpolation space of Tp'^ and Tp^ is Span{l, x, y. xy, y"^}. 

Define RQ = Span{l, X, y, 'z, xP' — y"^, xP' — z'^}. For the reference element K, let Rk 
be a interpolation operator i?g : H'^{K) —f RQ defined by 

%0(M,) = J(M,),j = l,...,6 

for all e i/^(K). It is exactly the so-called rotation element. Obviously, 



where 



i=l 



^1 = -gy -g^ 

■02 = — y z H — X a; H — , 

^ 6^ 6 3 2 6' 

^ U2 1^2^1^2^1^^1 

= -6^ +3^ +2^+6' 

= — x z H — y y^ — , 

^ 6 6 3^ 2^ 6' 

6 6 3 2 6 

= -6^ "6^ +r +2^+6- 



Thus we have 

6 



i=l 
6 



1=3 

+j(M2) (-^r - 



= e{K,<P,y,z) + <P{M,). 
It is easy to notice that ipi{l,y,z) = ipi{—l,y,'z) for i = 3, 4, 5, 6, which will leads to 

6 

1=1 

= ^ 0(M.)^^«(i, ^) + mi) (-^r - ^?') 

For any given K G we can define the interpolation operator i?K ^k'^k^- Denote 
by and the centroids of the faces with outward normal (1, 0, 0) and (—1, 0, 0), 
respectively. Then for any e if^(K), we have 

i?K0U=i = e(K, 0, z) + 0(M^+), 
i2K0|x=-i = e(K, 0, y, ^) + 0(M^-). 

3.3 The error estimates 

Let us give an estimation of and the estimation of Ey and E^ can be done 
similarly. According to the definition of E^ and Tp'^ ,Tp~ , we have 



Ex — ^ / v^otVujwds 



= ( / i/JaVuj(u) - T^+(u)))ds+ / i/JaVuj(«; - r^"H)ds 

Next we want to prove 

(«;|^.+ - n-'i^n)) = (^/;|^.- - H), (13) 
which can be done on the reference element. That is, we need to show 

- Tf{^)) = i^lpi- - Tfi^))- (14) 



Due to w e P^^, we have 



14 

w 

i=l 



Direct computation will derive 

(?.!]?!+ - Tf'm = (0^If!- - Tfi^^)), (15) 



which will ensure Eq. f ll4p . Denote by Mp{u) the value of u^cxWu at the centroid of 



3 

F. Thus we arrive at 



+ 



Ke. 
+ 

Kg 



J (u^aVu - Q{K,u,y, z) - Mp.^^{u)^ {w - T^-{w))ds^ 

J2 i I {i^'^cxVu - R^u) {w - T^-^{w))ds 

+ / {u^aVu- R-ku){w -T'^-{w))ds\ 
JFk- J 

where -Rkw is an interpolation polynomial of u'^cxWu. Here we remark that in the 
second step —Mp{u) is added due to the orthogonality ([8]) and in the third step 
— 0(K, M, ?/, 2;) is added due to Eq. (ITll) . We also mention that Rj^u can not be added 
directly according to Eq. (IT^ because Rj^u takes different value on the face Fp' and 
Since /?k and Tp preserves -Pi(K) and Pi(Ff ), respectively, it follows from trace 
theorem and Cauchy-Schwartz inequality, we get 

\Ex\ < Ch?\\w\\h\\u\\m(^n)- 

Similarly, we also have 

,2ii„..ii ii„.ii „ I zji I / ^r,2i 



Hence 



\Ey\ < Ch \\w\\h\\u\\H3{n), \E^\ < Ch \\w\\h\\u\\H3^n). 



\ah{u,w) - {f,w)\ _ \^£±3L±3i < Ch^Wul 



sup = sup j- S ^IT' \\U\\h3{Q)- 

>&.A"^^ IpIU wG.yrv^ IpIU 



By collecting the above results, we get the following energy-norm error estimate. 

Theorem 1. Let u e H^{^1) f] H^{^1) and Uh G satisfy ([3]) and (H]), respectively. 

Then we have the energy norm error estimate: 

\\u - Uh\\h < Ch'^\\u\\H3{n)- 
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By the standard Aubin-Nitsche duality argument, the L2(n)-error estimate can be 
easily obtained, but the corresponding proof is omitted. We state the result in the 
following theorem. 

Theorem 2. Let u e H^{Q) n i^o(^) "''^d Uh G be the solution of (E]) and 

respectively. Then we have 

\\u - Uh\\o < Ch^\\u\\H3{n)- 



4 Error estimates of other Smith-Kidger element 



In this section we claim that solving second-order elliptic problem with Smith-Kidger 
element of type 2 and 5 is also convergent with optimal order. In these case, it is easy 
to check that the orthogonality in Lemma |2] holds. The difference during the proof 
lies in the construction of the interpolation operator. For the second type element, 
the interpolation of T^, T^, should be Span{l, y, z, yz, y"^}, Span{l, x, z, xz, z^} and 
Span{l, xy, x^}, respectively. And for the fifth type element, the corresponding 
interpolation spaces should be taken as Span{l, y, z, yz, z^ + y^}, Span{l, x, z, xz, x^ + 
z^} and Span{l, x, y, xy, x^ + y^}, respectively. 

For Type 6 element, the orthogonality in Lemma |2] does not hold, but the Eq. ( IT4l) 
hods. Thus, we have 



I -Err 



+ [ {ly'^aVu- P^{iy'^aVu)){w-T^-{w))ds] 



< Ch\\u\\2\\w\\h 



where 



P°(i/^q:V'u) 



|K| 



ly aVuds 



K 



and |K| = jj^ds. By a similar derivation, we will get 

Theorem 3. Let u G H'^{Q) nifo(^) satisfy ([3]) and Uh he the solution of (jlj) with the 
sixth type element. Then we have the energy norm error estimate: 

\\u - UhWh < Ch\\u\\H2{n), 
\\u - Uh\\o < Ch'^\\u\\H2(^n)- 
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5 A new 14-node brick element 



In this section, we present a new element with 14-node. The degrees of freedom are 
the same with those in Smith-Kidger element and Meng-Sheen-Luo-Kim element. But 
the shape function space is taken as ® Span{x(y2 + z'^),y{x^ + z'^),z{x'^ + y^)}. The 
space can match the DOFs and the basis functions are 



en 



1 

16 



1 

16 



1 

16 



1 

16 



1 

16 



1 

16 



1 

16 



1 

16 



H — {xy + xz + yz + xyz) , 
8 

-1 + + + + x{y^ + z^) + y{x^ + z^) - z{x^ + y'^)) 

+^ {xy -xz-yz- xyz) , 
-1 + x'^ + y^ + z'^ + x{y^ + z^) - y{x^ + z^) + z{x^ + y^)) 

+- {—xy + xz — yz — xyz) , 
-l + x'^ + y^ + z'^ + x{y^ + z^) - y{x^ + z^) - z{x^ + y^)) 

+- {—xy — xz + yz + xyz) , 
-1 + x'^ + y^ + z^ - x{y^ + z^) + y{x^ + z^) + z{x^ + y^)) 

+- {—xy — xz + yz — xyz) , 
-1 + + + - + z^) + y{x'^ + z^) - z{x^ + y^)) 

+^ (~^2/ + - 1/2 + xyz) , 
-1 + + + 2^ - x(y^ + z'^) - y{x^ + z^) + z{x^ + y^)) 

+^ (xy - X2; - y^; + xy2;) , 
-1 + x^ + y^ + 2;^ - x(y^ + z^) - y(x^ + z^) - z{x^ + y^)) 

+^ (xy + xz + yz - xyz) , 



^12 = I - + 1 (-2:' + 2/' - + \y{x' + Z^), 
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^14 = ^ - + i {-x" -y^ + z^) + ^z{x^ + y^). 

To analyze the convergence, we need to verify the orthogonahty in Lemma [2] and Eq. 
f|T^ . The orthogonahty can be checked directly like the proof in Lemma [2l To satisfy 
Eq. f|T^ . we take the corresponding interpolation spaces as Span{l, + ?/^}, 

Span{l, X, z, xz, x"^ + z"^} and Span{l, x, y, x?/, + y^}, respectively. Thus, by the idea 
above, we can also get optimal convergence order for the second-order elliptic problems. 
That is, in this case, theorems [T] and |5] hold. 

6 Further remarks and conclusions 

In this paper, we have proved that for second-order elliptic problems, the Smith- 
Kidger element of type 1, 2 and 5 can obtain optimal convergence order both in energy 
norm and L2(fi) norm, while the sixth type element loses one order of accuracy in each 
norm. In the proof, the key points lie in that they have weak orthogonality (Lemma 
[2]) and satisfy Eq. f|T^ . In [5], we also proposed another kind of DOFs, that is, the 
value at the eight vertices and the integration value of six faces. Indeed, it is easy to 
check that Type 1, 2, 5 and the new element presented in this paper can get optimal 
convergence order for second-order elliptic problems. Besides, we can show that Type 6 
also has the optimal convergence order because the weak orthogonality naturally hods, 
which improve one order accuracy in this case. 
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